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Abstract. The composition of ice grains provides an important tool for the study of the molecular environment of star forming 
regions. Using ISAAC at the VLT to obtain spectra around 4.65 fim we have detected for the first time 'XCN' and CO ice in an 
extragalactic environment: the nuclear region of the nearby dusty starburst/ AGN galaxy NGC 4945. The profile of the solid CO 
band reveals the importance of thermal processing of the ice while the prominence of the XCN band attests to the importance 
of energetic processing of the ice by FUV radiation and/or energetic particles. In analogy to the processing of ices by embedded 
protostars in our Galaxy, we attribute the processing of the ices in the center of NGC 4945 to ongoing massive star formation. 
Our M-band spectrum also shows strong HIPfy6 and 0-0 S(9) line emission and gas phase CO absorption lines. The HI, H2, 
PAH, gas phase CO and the ices seem to be embedded in a rotating molecular disk which is undergoing vigorous star formation. 
Recently, strong OCN" absorption has been detected in the spectrum of the Galactic center star GC:IRS 19. The most likely 
environment for the OCN" absorption is the strongly UV-exposed GC molecular ring. The presence of processed ice in the 
center of NGC 4945 and our Galactic center leads us to believe that processed ice may be a common characteristic of dense 
molecular material in star forming galactic nuclei. 
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1. Introduction 

It has long been recognized that dust is an important com- 
ponent of the ISM in star forming regions in our galaxy and 
nearby galaxies. Optical studies of high redshift galaxies as 
well as deep mid-infrared and submm surveys have recently 
stressed the importance of dust in galaxies up to high red- 
shifts. Understanding the composition, origin and evolution 
of dust - particularly in star forming regions - is therefore a 
key question of astrophysics. Of special importance are inter- 
stellar ices, which are only present in molecular cloud envi- 
ronments. Mid-infrared spectra taken with the spectrometers 
aboard the Infrared Space Observatory (ISO) have shown ices 
to exist in a variety of extragalactic sources, from nearby star- 
burst nuclei up to distant dust-enshrouded ultra-luminous in- 
frared galaxies (ULIRGs). Weak H2O ice absorptions were first 
detected in the nuclear spectra of the nearby galaxies M 82 and 
NGC 253 (Sturm et al. 2000 1. Strong ice absorptions of H2O, 
CO2 and the XCN/CO blend were first seen in the nucleus of 
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NGC 4945 (Spoon et al. l2000b . This was followed by the dis- 
covery of H2O, 'HAC and CH4 ice in the nuclear spectrum of 
NGC 4418 (Spoon et al. 2001j. So far ices have been found in 
some twenty galaxies (Spoon et al. I2002I I. the most distant of 
which is IRAS 00183-71 11, a ULIRG at z=0.33 (Tran et al. 
.2001J . 

Interstellar ices are ideal probes for the conditions in the 
coldest and best shielded galaxy components - their molecu- 
lar clouds. Embedded protostars can process interstellar ice in 
their environment thermally as well as through FUV photoly- 
sis. Solid CO and 'XCN' provide prime probes for this pro- 
cessing. The fundamental vibrational modes of these species 
correspond to wavelengths around 4.65 jum in the M-band at- 
mospheric window. CO ice is highly sensitive to thermal pro- 
cessing of its environment. In quiescent molecular clouds, as 
probed by e.g. the field star Elias 16, the CO ice feature is 
dominated by the feature characteristic for an apolar ice mix- 
ture. When heated to temperatures above 20 K, this ice mix- 
ture sublimates and the only CO ice mixture surviving is a 
mixture dominated by H2O ice. This mixture has a distinctly 
different profile and sublimates at temperatures of ~90 K. The 
spectrum of the deeply embedded massive protostar W33A 
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is dominated by this type of CO profile (Chiar et al. I1998> . 
The presence of 'XCN', or OCN" ice, after its most popular 
identification (Demyk et al. 199S'), is indicative of strong pro- 
cessing of icy grains by UV light, energetic particles or heat- 
ing (Lacy et al H^S?! Grim et al]T^M Palumbo et al. 2000 
F. van Broekhuizen, priv. comm.). In our Galaxy, OCN" ice 
is strongest in the embedded massive protostar W33A (e.g. 
Pendleton et al. 1999). The feature is generally weak (com- 
pared to CO ice) in other protostars (e.g. Pendleton et al. 1999) 
and absent in quiescent molecular clouds, such as probed by the 
line of sight towards the field star Elias 16 (Chiar et al.[T995 1. 

Here we report on our VLT/ISAAC L and M-band follow- 
up spectroscopy of the rich ice absorption spectrum of the nu- 
cleus of NGC 4945 (Spoon et al. .2000) , a nearby (3.9 Mpc; 
l"=18pc; Bergman et al. \TW^ luminous (Lir=3x10'°Lo) in- 
frared galaxy, seen nearly edge-on (i~78°; Abies et al. 19871. 
The central region of this galaxy is dominated by a visually ob- 
scured starburst (Moorwood et al. 119961 Marconi et al. I2000t 
and a heavily enshrouded AGN, only seen in hard X-rays 
(Iwasawa et al. [T993l Guainazzi et al. l2000t . Paa (1.88yum) 
and K-band (2.2 urn) images obtained with HST/NICMOS 
(Marconi et al. I2000t reveal a complicated nuclear morphol- 
ogy, resulting from a nuclear starburst partially obscured by 
a strongly absorbing circumnuclear star forming ring, seen 
nearly edge-on. Fig.H shows the HST/NICMOS K-band im- 
age (Marconi et al. .2000j . rotated such that the galaxy major 
axis Ues horizontally. A possible geometry for the circumnu- 
clear starburst ring, seen under an inclination of -78°, is in- 
dicated by a dotted ckcle with radius 5.5" (lOOpc). The ap- 
parent asymmetric distribution of the nuclear emission interior 
to this ring (strong emission at position E, no emission at po- 
sition B) has been attributed to patchy absorption within the 
ring (Marconi et al. l2000> . Note that the emission at position A 
on the galaxy major axis appears to arise in the circumnuclear 
ring, not in the nuclear starburst. The dark structures seen at po- 
sition D in Fig.[nare likely gas filaments rising above the star 
forming ring and appearing in absorption against the bright nu- 
clear continuum. The Paa image (Marconi et al. 2000J shows 
more of these dark filamentary structures, one of which passes 
in front of the K-band nucleus. The AGN, the second bright- 
est Seyfert 2 nucleus in the sky at hard X-rays (Done et al. 
I1996> . does not appear in any of the HST/NICMOS images. 
Groundbased lOyum ESO/TIMMI images (1" seeing) reveal 
no sign of the central monster either (Marconi et al. 2000 1. 
Also mid-infrared spectroscopic AGN tracers, like 7.65 /im 
[Nevi] and 14.3 & 24.3 //m [Nev], resulted in strong upper 
Umits (Spoon et al. .2000) . It is hence Ukely that the AGN is 
strongly obscured in all directions by material which is most 
probably close to the black hole for the obscuration to be ef- 
fective. The conical cavity, protruding from the nucleus along 
the minor axis and seen in the near-infrared (Moorwood et al. 
[19961 and in soft X -rays (Schurch et al. 2002 1, is therefore not 
a Seyfert ionization cone, but a cavity cleared by a supernova- 
driven starburst 'superwind' (Moorwood et al. ll996> . 

Here we present infrared L and M-band spectra of all nu- 
clear components discussed above obtained using VLT/ISAAC 
with its slit aligned with the galaxy major axis (Fig.[0. 



2. Observations 

The 3-5 fim spectra of the NGC 4945 nucleus were ob- 
tained using the Infrared Spectrometer And Array Camera 
(ISAAC) mounted on the VLT (UTl-Antu) at the ESO Paranal 
Observatory. The observations were performed on the night of 
June 15, 2000 for the L-band spectrum and on the nights of 
June 15-16, 2000 and July 1-3, 2001 for the M-band spec- 
trum. The low resolution grating and the 1" slit were used for 
the L-band spectrum with a corresponding resolving power of 
R=360. The M-band spectrum was obtained using the medium 
resolution grating and the 1" slit, resulting in a spectral resolv- 
ing power of R=3000 (Av=100km/s FWHM). In both cases, a 
single setting was done covering the entire L-band and the re- 
gion from 4.56 to 4.80//m in the M-band. The spectra have a 
total integration time of 240 s and 5123 s for the L and the M- 
band, respectively. The telescope was operated in the standard 
chop-nod mode with a chop throw of 20". The telluric standard 
BS 5571 (B2III) was observed just before or after each obser- 
vation of NGC 4945 with a maximal airmass difference of 0. 15. 

The data were reduced using our own IDL routines. The 
individual frames were corrected for the non-linearity of the 
Aladdin detector, distortion corrected using a star trace map 
and bad pixels and cosmic ray hits were removed before co- 
adding. For the M-band this procedure resulted in six co-added 
source frames and six co-added standard star frames. Standard 
star spectra were obtained by extracting the positive spec- 
tral trace from each of the six standard star frames. The six 
source frames were then divided by their associated standard 
star spectra, taking into account an optimal small shift between 
the source frame and standard spectrum by requiring that the 
pixel-to-pixel noise on the continuum of the final source frame 
be minimized. No correction for airmass differences was at- 
tempted due to insufficient signal-to-noise ratio of the source 
frames. No hydrogen recombination lines were detected in the 
standard spectrum and thus no attempts were made to correct 
the standard for photospheric lines. Next, each source frame 
was flux calibrated relative to the standard and wavelength cal- 
ibrated using the telluric absorption lines in the standard star 
spectrum. The final source frame was then produced by stack- 
ing the six source frames, taking into account small dispersion 
shifts among the frames. 

The flux calibration is estimated to be better than 15% and 
the wavelength calibration is accurate to 150km/s and 15 km/s 
for the L and M-band spectra, respectively. 

3. Results 

3.1. The Sfim water ice band 

The ISAAC L-band spectrum of the central region of 
NGC 4945 (Fig.EJ is dominated by a broad (2.1-4.0 fim) ab- 
sorption band attributed mainly to water ice (e.g. Smith et 
al. [TWyi Chiar et al. UUU^ . The width of the feature is best 
appreciated in the low resolution 24"x24" ISO-PHT-S spec- 
trum, which is shown for comparison. Superimposed on the 
absorbed L-band continuum are the PAH emission bands at 
3.3 & 3.4jum as well as the H Bra line at 4.05 jt/m. The ISAAC 
spectra of the central 2"xl" and 11.5"xr' agree weU with 



H.W.W. Spoon et al.: Detection of strongly processed ice in the central starburst of NGC 4945 



3 




10 5 0-5 -10 

Offset [arcsec] 

Fig. 1. HST/NICMOS K-band image of the nucleus of NGC 4945 (reproduced from Marconi et al. (2000 1. The image has been 
rotated so that the galaxy major axis (PA=43°) is oriented horizontally. The position of the VLT/ISAAC slit is marked by crosses, 
spaced 1" apart. Arrows indicate several positions of interest, which are discussed in the text. The dotted ring gives an impression 
of how a circumnuclear ring with radius 5.5" (100 pc) would be seen under an inclination of 78°. 



the larger beam ISO-PHT-S spectrum, assuming the water 
ice feature in the latter spectrum to be diluted by a non- 
absorbed stellar and dust continuum. Adopting a flat contin- 
uum (fixed at 4.0^. 1 yum) and depending on whether we fit 
the bottom or the blue wing of the feature, we find a water 
ice column of 41-47 x 10'^ cm"^ (assuming a band strength of 
2.0x10"' cm/molecule; Gerakines et al. 1995 1. Note the ab- 
sence of 3.94jum [Siix] line emission from the ISAAC L- 
band spectrum. This line arises in soft X-ray photoionized gas 
and is comparable to or brighter than the H Bra line in many 
Seyfert galaxies (Oliva et al. [T994l Lutz et al. l2002l l. The ab- 
sence ([Si ix]/HBrff <0.10) may be taken as evidence for the 
extremely high obscuration of the AGN in NGC 4945 in the 
L-band. 

3.2. Processed CO and OCN- ice 

The M-band spectrum of the central 2"xl" of the nuclear re- 
gion of NGC 4945 (Fig.|3} contains strong and relatively broad 
absorption features at 4.62 and 4.67 fim as well as gas lines due 
to H, H2 and CO. The line centers of the gas phase lines shift 
as a function of position along the slit. 

In order to investigate the various M-band absorption fea- 
tures, we fitted a linear continuum to the two pivot ranges 4.54- 
4.55 yum and 4. 78-4. 79 /^m. The resulting optical depth spec- 
trum is shown in Fig.|4t. In Fig.^J),^ and^ we show spectra 
of three comparison objects: the embedded massive protostar 
W33A with the strongest known XCN feature; the Galactic 
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Fig. 2. The VLT/ISAAC rest frame L-band spectrum of the 
central 2"xl" (36xl8pc2) and 11.5"xl" (210xl8pc2) of 
NGC 4945 (black) compared with the ISO-PHT-S spectrum 
obtained in a 24"x24" aperture (grey histogram). The arrow 
indicates the expected wavelength of the 3.93 yum [Si ix] line. 



Center (Moneti et al. 1200 1 > showing XCN absorption towards 
its massive star forming region which is obscured in what is 
eff'ectively an 'edge-on' view towards the center of our own 
galaxy; and, finally, the unprocessed line of sight to the field 
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Fig. 3. The VLT/ISAAC rest frame M-band spectrum of the 
central 2" x 1 " (36 x 1 8 pc^) of NGC 4945 . The spectrum shows 
emission lines of HIPfjS and H2 0-0 S(9), broad absorption 
features of ices containing XCN and CO, and absorption lines 
of gas phase ^^CO and '^CO. The adopted continuum and the 
shape of the ice feature, corrected for the presence of gas phase 
'^CO and '^CO lines, are drawn in grey. 



star Elias 16 located behind the Taurus molecular cloud (Chiar 
et al.[T995}. The similarity of NGC 4945 to W 33 A is particu- 
larly striking and will be discussed later. 



Following Galactic ice studies (Tielens et al. ll991l Schutte 



& Greenberg 119971 Ch iar et al. 119981 Demyk et al. 119981 
Pendleton et al. '1999'i we identify the prominent absorp- 
tion feature centered at 2168.0cm"' (FWHM=24.0cm"') with 
soUd state 'XCN', or OCN" after its most likely identification 
(Demyk et al. '19981 In order to measure the shape and depth of 
the feature, we adopt the ice feature shown in Fig.|3] This con- 
tinuum has been corrected for the presence of gas phase '^CO 
and '^CO absorption lines, discussed later in this Section. The 
OCN" feature appears slightly blueshifted and narrower than in 
W 33A (vo=2165.5 cm"' ; FWHM=26.7 cm"'). The differences 
are similar to those found between the embedded protostars 
AFGL961 and W 33 A (Pendleton et al. ,1999 l and may be at- 
tributed to differences in ice composition (Grim & Greenberg 
[T987t . The column density of OCN" is Hsted in Table[T]and 
were computed assuming an OCN" bandstrength of 1.3x10"' 
(F. van Broekhuizen, priv. comm.). 

Detailed studies of the observed profiles of the solid 
CO ice band in Galactic sources have shown that they can 
all be fitted by varying contributions of the same three 
components (Boogert et al. I2002all2'002 b" Pontoppidan et 
al. .2003J . These three components are thought to repre- 
sent "pure" solid CO (Gaussian centered at V()=2139.9cm ' 
with FWHM=3.5cm '), traces of CO trapped in H2O ice 
("polar CO"; Lorentzian centered at vo=2136.5cm ' with 
FWHM=10.6cm-i) and CO in soUd CO2 ice ("apolar CO"; 
Gaussian centered at vo=2143.7cm"' with FWHM=3.0cm"') 
(Boogert et al. i2002at Pontoppidan et al. .2003J . Here we note 




4.55 



4.75 



4.60 4.65 4.70 

Rest wavelength [fj.m] 

Fig. 4. The optical depth spectrum of ices in NGC 4945 
compared to ices seen in Galactic lines of sight. The Top 
panel shows the optical depth spectrum of the central 2"xl" 
(36xl8pc2) of NGC 4945 (grey), determined adopting the 
continuum shown in Fig.|3] Also indicated are the results of the 
fit to the ice profile. The OCN" fit component is denoted by a 
dashed line, the polar CO component by a dotted line, the pure 
CO component by a continuous line and the apolar compo- 
nent by a grey surface. The overall fit is shown as a black line. 
The spectrum was obtained at a resolving power of R=3000. 
Second panel: Optical depth spectrum for the massive embed- 
ded protostar W 33A. The spectrum has been degraded to a re- 
solving power of R=750. Third panel: Idem for the ISO-SWS 
line of sight to the Galactic Center (Sgr A*). The spectrum was 
obtained at a resolving power of R=2000. Bottom panel: Idem 
for the unprocessed line of sight to field star Elias 16, located 
behind the Taurus molecular cloud. The spectrum was obtained 
at a resolving power of R=1200. 



that the spectrum of NGC 4945 is dominated by traces of CO 
trapped in H2O ice with little or no evidence for the other 
two components (Fig.0t). In contrast, the spectrum of W 33A 
(Fig.|4j)) has a noticable contribution of the "pure" solid CO 
component (Tielens et al. 1991'!. Table^lists the column den- 
sities of the various CO ice components, which were com- 
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puted assuming a CO bandstrength of 1.1 xlO'^ cm/molecule 
(Gerakines et al. 1995). 

The nuclear spectrum of NGC 4945 also reveals CO in the 
gas phase. In Fig.|3]we identify a total of 13 fundamental ro- 
vibrational absorption lines of '^CO and several of '^CO. This 
is a sufficient number of lines to attempt a single component 
model fit to determine the temperature, intrinsic line width and 
'^CO and '^CO gas column densities. For this purpose we use 
the isothermal plane-parallel LTE CO gas models of Cami et 
al. (2002), folded with the appropriate VLT/ISAAC spectral 
resolution (R=3000; Av=100km/s FWHM). In the fitting pro- 
cedure both the observed spectrum and the model spectra are 
normalized, through division by a strongly smoothed version 
of the respective spectra. A four-parameter least-squares min- 
imalization procedure then picks the best fitting model. Fig.|5l 
shows the result for the nuclear spectrum of NGC 4945. The 
CO gas appears moderately warm at 35^2 5 with an intrin- 
sic line width (FWHM) of 50+5 km/s and column densities log 
N(i2cO)=18.3+0.1cm-2 and log N(i3CO)=17.2+[J (^^ ^,^-2^ 

suming a covering factor 1 for the absorber. The uncertainties 
listed above do not take into account the systematic errors re- 
sulting from the data reduction, which may be appreciable. The 
very low '^CO/'^CO ratio of ~13 indicates that some '^CO 
lines are likely to be optically thick. In addition, the absorb- 
ing material may not be in a uniform screen 'covering fac- 
tor r configuration (Fig.^. This will particularly affect the 
derived column for '^CO. A more realistic, yet more uncer- 
tain, value for the '^CO gas column may be derived from the 
'^CO column. Assuming a '^CO/'-'CO ratio of 80, we find log 
N('^CO)=19.1+gg5 cm"^. Note that the (high resolution) spec- 
trum of the physically closest resembling object, the protostar 
W 33A, also contains CO gas phase Unes (Mitchell et al. 1988 ). 
At the resolution of the ISO-SWS spectrum (R~750; Gibb et 
al. 1^000) these are however not detectable. 



3.3. The spatial distribution of tlie ice 

The analysis presented above focusses on the nuclear spectrum 
of NGC 4945, which comprises the central 2"xl". The nuclear 
continuum emission extends however from ~4.5" SW to ~2" 
NE of the nucleus, measured in a 1" wide slit along the SW-NE 
oriented (PA=43°) galaxy major axis (see Fig.Q. Line emis- 
sion (Fig.0 can be traced as far out as -4.5" SW (position E) 
and ~5.5" NE (position A) of the nucleus. Both the continuum 
and line emission appear clearly weakened in a ~1" long strip 
located 1.7" SW (position D) of the nucleus. This strip coin- 
cides with a dark filamentary structure in the HST/NICMOS K- 
band image (Fig.Q of Marconi et al. (I2OOOI ). The same image 
further shows strong extinction from what might be a circum- 
nuclear ring to be responsible for limiting the observable part of 
the major axis continuum to the range found (see Sect. 1). Our 
analysis further shows that the ice and gas absorption features 
are detected wherever there is background continuum to ab- 
sorb. The depth of the ice features appears to be quite constant 
(Fig.|6j, except for position D, where the solid state absorptions 
may be significantly stronger. Again interestingly, this posi- 
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Fig. 5. CO gas absorption lines in the spectrum of the central 
2"xl" (36 X 18 pc^) of NGC 4945. The observed spectrum, af- 
ter dividing out the ice feature, is shown in black. The best 
fitting CO gas model is drawn in dark grey. Fit residuals are 
shown in light grey. 



tion coincides with the dark filamentary structure seen in the 
HST/NICMOS K-band image (Fig.Q. 

The lower panels of Fig.0show position-velocity diagrams 
for the HI PfjS and H2 0-0 S(9) emission Unes. The distributions 
are remarkably similar over the central +2", both indicating ro- 
tation about the nucleus. The only significant deviation occurs 
at a position 0.7 " NE of the nucleus (position C), where the H2 
0-0 S(9) emission extends -40 km/s beyond the highest veloc- 
ity traced by HI Pf/3. We speculate that this position may coin- 
cide with the start of the N-S oriented molecular ridge, traced 
in H2 1-0 S(l) by Moorwood et al. (I1996> and Marconi et al. 
(I2OOOJ . bordering the conical cavity cleared by a starburst su- 
perwind (Moorwood et al. lI996> . Another deviation occurs 3"- 
4" NE of the nucleus (position B), where a dark cloud in the 
circumnuclear ring weakens the line emission from the nuclear 
starburst it eclipses. Given the high inclination of the circum- 
nuclear ring, the dark cloud may actually be in an orbit in the 
outer part of the circumnuclear ring, with an appreciable tan- 
gential orbital velocity component. This would explain why the 
HI Pf/3 emission from that cloud is ~50 km/s lower than that of 
the bright patch 1 " NE of it (position A). The top panel of Fig.0 
shows the velocity curves derived from fitting single Gaussians 
to the line emission in each of the observed HI Pf/3 and H2 0- 
S(9) position- velocity distributions. The velocity curves ob- 
tained for both emission lines are consistent and yield the same 
velocity gradient of ~60km/s/arcsec across the nuclear region. 
The heliocentric systemic velocity of 56 1+4 km/s (Dahlem et 
al. ll993> occurs at the brightest M-band cross dispersion pixel. 
Away from the nucleus the velocity curves flatten ofi" to rota- 
tional velocities 120±10km/s above and below systemic. This 
value agrees well with the results for HBry and H2 1-0 S(l) 
obtained by Moorwood & Oliva (1994'). Note that Ott et al. 
(.2001 J find rotational velocities up to 160 km/s for HI (21cm). 
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Fig. 6. Optical depth spectra of ice features in NGC 4945 for 5 adjacent positions (3.2" SW: triple dot-dashed, 1 .7" SW: dotted, 
0.45" SW: solid, 0.60" NE: dot-dashed, 1.6" NE: dashed) along the galaxy major axis. The spectra have been smoothed and 
rebinned to lower spectral resolution. The optical depth spectra reveal no significant differences, except for the generally larger 
optical depth at position D, 1.7" SW of the nucleus (dotted line). Left: Optical depth spectra of water ice. The feature centered 
at 3.3 jum is due to the presence of PAH emission in the water ice absorption feature. Right: Optical depth spectra of XCN and 
CO ice. 



Assuming that the measured velocities trace the potential, the 
mass within the nuclear region amounts to 1x10** Mq at 50 pc 
and 3 xIO'^Mq at lOOpc. Also shown in the upper panel of 
Fig.0is the velocity curve for gas phase CO. The CO absorp- 
tion line velocities were derived by cross-correlating the CO 
absorption line spectra with the best fitting CO gas model. In 
order to acquire sufficient signal-to-noise in the lines, the spa- 
tial information has been combined into 5 bins of 1.0"- 1.5" 
along the major axis. The resulting velocity curve for the CO 
absorption lines (excluding the point at 1.7" SW; position D; 
see below) clearly samples a different velocity field than the 
emission lines do. Based on the four CO points, we measure a 
velocity gradient for the CO gas of 17 km/s/arcsec across the 
nucleus. The velocity measured for the 1.7" SW point (posi- 
tion D) clearly does not fit in with the velocity gradient seen 
for the other four CO points. Interestingly, the HST/NICMOS 
K-band image (Fig.Q shows at this position a dark filamen- 
tary structure against the bright K-band continuum. Possibly, 
the CO gas radial velocity measured at 1.7" SW is related to 
this foreground (circumnuclear ring?) structure instead of to 
the rotating inner ring/disk. 

4. Discussion 

The wealth of ISM features detected in the VLT/ISAAC 
L & M-band spectra of NGC 4945 allows us for the first time to 
study simultaneously the conditions of ionized hydrogen gas, 
molecular hydrogen, PAHs, icy grains and cold CO gas in the 
central region of another galaxy. Combined with the spatial in- 
formation obtained in a I" wide slit oriented along the galaxy 
major axis, our observations give insight in the kinematics and 
spatial location of the different ISM components probed. 

The picture that emerges is of an extended, fragmented and 
vigorously star forming rotating molecular ring/disk, surround- 
ing the deeply buried, and hence passified, AGN. The H, H2 and 
PAHs trace the interaction of the massive stars - interspersed 



Table 1. Measured column densities (10'^ molecules/cm^) and 
column density ratios for NGC 4945, the massive embedded 
protostar W 33A, the line of sight towards the Galactic center 
(Sgr A*) and towards a field star, Elias 16, located behind the 
Taurus molecular cloud. 





N4945 


W33A 


Sgr A* 


Elias 16 


N(polar CO ice) 


9.4 


11.4 


2.8' 


2.9 


N(pure CO ice) 


0.19 


1.4 




4.3 


N(apolar CO ice) 


0.16 


0.0 




1.0 


N(total CO ice) 


9.7 


12.7 


3.5" 


8.2 


N(XCN ice) 


1.6 


2.6 


0.2" 


<0.13'' 


N(cold CO gas) " 




130" 


66'' 


18'' 


N(warm CO gas) " 


130 


140" 


Id 




N(CO gas) 


130 


270" 


6T' 


18'' 


N(H20 ice) 


41-47 


110* 


llf 


25« 


N(C02 ice) 


>2.0' 


14^ 


IJj 


5^ 


N(polar CO)/N(CO) 


0.97 


0.89 


0.8 


0.35 


N(XCN)/N(CO) 


0.17 


0.20 


0.07 


<0.02 


N(C0)/N(H20) 


0.21-0.24 


0.12 


0.3 


0.33 


N(XCN)/N(H20) 


0.034-0.039 


0.024 


0.02 


<0.005 


N(CO ice)/N(CO gas) 


0.08 


0.047 


0.05 


0.5 



a) We define CO gas as cold if T„as <30 K and as warm if Tgas >30 K; 

b) Gibb et al. ROT : c) Mitchellet al. 1 198|(; d) Moneti et al. nMT\ : 
e) Based on the Sgr A* spectrum of Moneti et al. <2001> : f) Chiar et 
al. ( 2000 (; g) Chiar et al. (19951 : h) Whittet et al. iT989l : i) Spoon et 
al. CTUt obtained in a 24"x24" aperture; J) Gerakines et al. 119991 



within the ring, as indicated by the velocity behavior of the H 
and H2 lines - creating HII regions and PDRs. The extended 
continuum emission source created this way, forms the back- 
ground against which we see the various absorption features. 
The contribution of the AGN to this continuum is insignificant, 
as may be concluded from the absence of a pointlike source 
in both the HST K-band image (Marconi et al. ^000 1 and our 
ISAAC spectra. Hard X-ray observations have shown that the 
direct view of the AGN is blocked by a hydrogen column den- 
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sity of lO^-*^ (Iwasawa et al. [T^ Ay -2500). The ab- 
sence of the 3.94//m [Siix] coronal line and other NLR fine 
structure lines (Fig.|2 Spoon et al. 2000) from this region in- 
dicates that the ionizing radiation does not escape in other di- 
rections either. The shallow velocity curve of the cold CO gas 
is consistent with the CO being distributed in the outskirts of 
the nearly edge-on, rotating molecular ring/disk, of which the 
CO gas velocities are sampled only in front of the bright nu- 
clear continuum, where the line of sight velocity component 
increases only slowly with projected distance (see Fig.|8|i. An 
indication of its scale can be obtained in the simplified assump- 
tion of absorption in a single ring rotating at 120km/s. To re- 
produce the observed velocity gradient of CO absorption, the 
radius of this ring would be ~7" (~120pc). Other distributions 
of absorbing material on similar scales, and in particular similar 
to the dust ring seen in the data of Marconi et al. (2000 ), will re- 
produce the data as well. The material seen in HI absorption by 
Ott et al. ( 1200 U and the rotating central (<20") mm-wave CO 
emission likely belong to the same structure. Although there is 
no velocity information, we are tempted to co-locate the H2O, 
CO and OCN^ ice with the CO gas. First, the solid H2O, CO 
and OCN" can be traced as far out as the gas phase CO (Fig.|6j. 
Second, the optical depth of the solid state features is constant 
over the spatial range probed (except for the 1.7" SW point in 
Fig-Ell, implying co-location of the ices in a foreground posi- 
tion. Third, the filamentary structure seen at 1.7" SW (posi- 
tion D) shows both deeper solid state features (Fig.|6} and a de- 
viant CO gas velocity, suggesting co-location of the ices with 
the CO gas. Fourth, the absence of the pure CO ice compo- 
nent indicates grains which have been thermally processed to 
above 20 K. Hence, unlike for dark cloud lines of sight in our 
galaxy (such as Elias 16 in the Taurus molecular cloud), these 
ice grains are not located in some random, foreground, dark 
cloud along the line of sight, but instead close to a source of 
thermal heating. Fifth, the deep OCN" ice band shows that the 
ice is heavily processed. That again locates the ices within re- 
gions of massive star formation or close to the AGN rather than 
some foreground material. From the spatial extent and unifor- 
mity of the OCN" ice absorption (Fig.|6|l, a location close to 
the AGN can, however, be excluded. The only likely location 
therefore remain within or close to regions of massive star for- 
mation in the circumnuclear starburst. Indeed, within our own 
galaxy, the massive protostar W 33A is an extreme example of 
processed ices (Chiar et al. lT^^ Gibb et al. HM^ . Dark cloud 
material in general does not show any evidence for the OCN" 
absorption band (Fig.|5l; Whittet et al. .2001j . In fact, no other 
Galactic, luminous protostar shows such a strong OCN" band. 

The presence of OCN" in ice grain mantles is often taken 
as a sign of energetic processing by particles or UV photons. 
Recent experiments suggest however that also thermal process- 
ing may result in the formation of OCN" (F. van Broekhuizen, 
priv. comm.). A fouth possibity presents itself in the vicinity 
of an AGN: processing by AGN X-ray photons. The AGN 
in NGC 4945 is, however, strongly obscured and only hard 
X-ray photons manage to escape. Assuming similar obscura- 
tion towards the r~100pc ice region in NGC 4945 as towards 
our line of sight, we estimate a hard X-ray (20-100 keV) flux 
of ~0.5 erg/cm^/s from the observations of Guainazzi et al. 



( 12000 ). This is two orders of magnitude less than the UV flux in 
a photon-dominated region next to massive stars (e.g. Tielens & 
Hollenbach ll985t . The available very hard (unabsorbed) AGN 
X-rays are hence energetically insignificant compared to the 
starburst UV photons. In addition, such hard X-rays may not 
couple efficiently into individual ice grains. Further laboratory 
studies are required to settle these issues. 

The starburst in NGC 4945 is similar in luminosity to 
the prototypical starbursts M82 and NGC 253 which are also 
located in the nuclear regions of almost edge-on galaxies. 
Comparing ISO spectra of NGC 4945 (Spoon et al. 20001) to 
ISO data for M 82 and NGC 253 (Sturm et al. 2000; Forster 
Schreiber et al. 2003.) . the obscuration of NGC 4945 is clearly 
higher towards the ionized medium (as derived from the 
18. 71/33.48 /vm [S ill] ratio in the low density limit), towards 
the PAH emitting medium (8.6 fim and 11.3 yum PAH emission 
features strongly suppressed by 9.7 //m silicate absorption), and 
in 3 jum water ice absorption (rice ~2.5). It remains unclear, 
however, whether this is just an on average higher absolute ab- 
sorbing column towards the nuclear region of NGC 4945, as 
perhaps not implausible given also the slightly smaller physical 
size of its starburst, or whether there are differences in the prop- 
erties of the absorbing medium. Higher quality M band spectra 
of M82 and NGC 253 are needed to test whether those star- 
bursts also host processed ices as NGC 4945, but with columns 
corresponding to their lower obscuration. 

The spectrum of the nearest edge-on galactic nucleus, our 
Galactic center, also shows absorption features due to CO and 
OCN^ ice. The spectrum of GCTRS 19 (Chiar et al. l2002t 
shows the features at similar relative strengths as in W 33A 
and NGC 4945. The star itself (M star) is most likely not the 
source of the energetic processing. A chance projection with 
the processed surroundings of a young foreground star is pos- 
sible, but also unlikely. The projected distance of GC:IRS 19 
to the GC molecular ring is, however, small. The material in 
this ring is exposed to high UV fluxes from the central cluster, 
crudely comparable to the radiation field within the molecular 
ring in NGC 4945. It is hence plausible that energetic process- 
ing has created similar ice properties in both NGC 4945 and 
in the GC molecular ring. The 'pencil beam' line of sight to 
GC:IRS 19 may pass through this processed ring material, ex- 
plaining the presence of strong OCN" absorption in the spec- 
trum of GC:IRS 19. The larger beam 14"x20" ISO-SWS spec- 
trum of Sgr A* also shows CO and OCN" ice (Fig.|2 Moneti 
et al. 2001 ). The OCN" column is, however, nearly an order of 
magnitude smaller than in the pencil beam towards GC:IRS 19. 
This difference is likely due to the large number of background 
stars within the ISO-SWS beam, whose combined line of sight 
results in the observed low OCN" optical depth. While some 
of these pencil beams pass through the processed GC molec- 
ular ring, other pencil beams only sample unprocessed fore- 
ground material. If the latter line of sight dominates within 
the 14"x20" ISO-SWS beam, this would be a natural expla- 
nation for the low OCN" optical depth in the Sgr A* spectrum 
of Moneti et al. (.2001.) . 

The presence of processed ice in the centers of NGC 4945 
and the Galaxy leads us to believe that processed ices are one 
characteristic of dense molecular material in star forming nu- 
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clear regions of galaxies, and can be detected in favorable ori- 
entations. The amount of processing may then be a measure of 
the nuclear star formation activity and/or geometry. 

5. Conclusions 

Using ISAAC at the VLT to obtain spectra at 2.85^.lOiJ.m 
and 4.55^.80/zm, we have detected for the first time 'XCN' 
and CO ice in an extragalactic environment: the central region 
of the nearby dusty starburst/ AGN galaxy NGC 4945. 

The profile of the solid CO band reveals the importance 
of thermal processing of the ice while the prominence of the 
XCN band attests to the importance of energetic processing of 
the ice by FUV radiation, energetic particles and/or heating. In 
analogy to the processing of ices by embedded protostars in our 
Galaxy, we attribute the processing of the ices in the center of 
NGC 4945 to ongoing massive star formation. 

Our M-band spectrum also shows strong HI Pf/3 and H2 0- 
S(9) line emission and gas phase CO absorption lines. The 
HI, H2, gas phase CO and the ices seem to be embedded in a 
rotating molecular disk which is undergoing vigorous star for- 
mation. 

The non-detection of the 3.94//m [Siix] coronal line in 
our VLT/ISAAC L-band spectrum is in full agreement with 
the very high obscuration towards the AGN derived from mid- 
infrared coronal line observations and from the absence of 
a point-like source in K-band and N-band images. With the 
source of the obscuration probably close to the AGN for the 
obscuration to be eff'ective, the radiation from the AGN cannot 
be responsible for processing the ices. 

The obscuration towards the starburst in NGC 4945 is far 
higher than towards similar starbursts in M 82 and NGC 253, 
also seen nearly edge-on. It is unclear whether this is just an on 
average higher absolute absorbing column towards the nuclear 
region of NGC 4945, or whether there are diff'erences in the 
properties of the absorbing medium. Higher quality M-band 
observations of M 82 and NGC 253 are required. 

Recently, strong OCN" absorption has been detected in 
the spectrum of the Galactic center star GC:IRS 19. The most 
likely environment for the OCN" absorption is the strongly 
UV-exposed GC molecular ring. The presence of processed ice 
in the center of NGC 4945 and our Galactic center leads us to 
believe that processed ice may be a common characteristic of 
dense molecular material in star forming galactic nuclei. 
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Fig. 7. Position-velocity information for several species de- 
tected in the M-band spectrum of NGC 4945. Top panel: 
Velocity curves for HIPf/3 (black) and H2 0-0 S(9) (grey) as 
determined from the position-velocity diagrams in the panels 
below. The data have been rebinned to half the size of a seeing 
element (0.5"). The five large dots denote the CO gas line ve- 
locities in five positional intervals along the galaxy major axis. 
A linear fit to the CO gas line points (excluding the point at 
1.7" SW, which seems to trace a separate filament) is shown 
as a dashed line. The position of the nucleus is indicated by 
the vertical dotted line, the systemic velocity (561 km/s) by the 
horizontal line. Middle panel: Position-velocity diagram along 
the major axis for the HI PfyS emission line. Contours run from 
4 to 1 Oct in steps of 2cr. The black points indicate the center 
of the gauss fit to the line profile at indicated positions along 
the galaxy major axis. The spatial and spectral resolution is in- 
dicated in the lower left corner. Lower panel: Idem for the H2 
0-0 S(9) emission hne. Contours run from 4 to 14 cr in steps of 

2cr. 
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Fig. 8. Schematic face-on view of the central region of 
NGC 4945. In our model the HIPf/3 and H2 0-0S(9) emis- 
sion lines trace the rotation of the inner disk. Their velocity 
curves (Fig.^t) indicate that the rotational velocity increases 
outward to become constant (~120km/s) at r~60pc. Assuming 
the CO gas phase absorption lines to arise in the outskirts 
of the rotating disk (Vrot=120km/s), the CO velocity gradient 
(17 km/s/arcsec ^1 km/s/pc; Fig.^) sampled in front of the nu- 
clear continuum indicates that the "CO ring" must have a radius 
of ~120pc. The white strip in the CO ring indicates the section 
of the CO ring where CO gas absorption lines can be detected 
against the nuclear continuum. The line of sight sight velocity 
component (Vrad) in this strip changes only slowly with pro- 
jected distance, in agreement with the observed shallow CO 
velocity curve (Fig.Et)- 



